Statics and dynamics of a harmonic oscillator 
coupled to a one-dimensional Ising system 
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Abstract. We investigate an oscillator linearly coupled with a one-dimensional Ising system. The 
coupling gives rise to drastic changes both in the oscillator statics and dynamics. Firstly, there 
appears a second order phase transition, with the oscillator stable rest position as its order parameter. 
Secondly, for fast spins, the oscillator dynamics is described by an effective equation with a 
nonlinear friction term that drives the oscillator towards the stable equilibrium state. 
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We analise the behaviour of a classical harmonic oscillator (mass m, frequency C0q, 
position x and momentum p) linearly coupled to a one dimensional Ising system (N 
spin variables o t • = ±1) in contact with a heat bath at temperature T. After a suitable 
nondimensionalization of the variables, the system energy per spin site is 
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je(x,p, a) = P —^- - 1 £ a^+i , (1) 

The oscillator dynamics follows from Hamilton's equations, x+x = (1/N) Y!u=i Oity+l, 
while the spins dynamics is governed by a master equation with Glauber single-flip 
transition rates D20. Then, the rate at which the z-th spin flips in configuration O is 
Wi(a\x,p) = 5 [1 - YCj{Cj-\ + o i+ i)/2] /2, with y = tanh(2x/0). In these rates, 5 gives 
the (dimensionless) characteristic attempt rate for the spin flips and is the dimen- 
sionless temperature [HI]. The equilibrium properties of the system are obtained from 
the canonical distribution 0P eq {x,p, a) exp [— NJ$?(x,p, g)/Q\. By summing over the 
spin variables, we calculate the marginal probability for the oscillator variables, 
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in which the spins produce a non-harmonic additional contribution to the effective 
potential of the oscillator, %ff(x). The equilibrium stable rest points x eq of the oscillator 
are the minima of so that they verify the bifurcation equation x eq — tanh(x e9 / 0) = 
0. Thus, there is a second order phase transition at 9 = 1 , with x eq as its order parameter. 

In the limit of fast spins, a perturbative analysis in powers of 8 1 <ti 1 yields an 
approximate nonlinear equation for the macroscopic (completely neglecting the fluctua- 




FIGURE 1. Averaged stochastic trajectories (x(t)) = x(t) (circles) versus nonlinear (solid blue line) 
predictions for (a) 9 = 4, (b) 9 = 0.95, and (c) 9 = 0.6. 
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The spins are so fast that they relax, almost instantaneously, to the equilibrium state 
corresponding to x(t) . But, since x(t) is slowly varying, there is a small separation from 
equilibrium proportional to dx/dt, which gives rise to the friction term. 

In the figure, the theoretical result © (its numerical integration) is compared to the 
simulation of the stochastic model, for a random initial configuration of the spins and 
suitable initial conditions for the oscillator [1]. For high temperatures, case (a), and for 
a temperature below (but close to) the critical value = 1, case (b), the agreement 
between simulation and theory is excellent. On the other hand, for a further lower value 
of the temperature, case (c), the theory breaks down. This was to be expected, since 
for very low temperatures the spins relaxation time diverge and they are actually slow 
as compared with the oscillator. Hence, another approach is needed, in which the fast 
vibrations of the oscillator are averaged using the method of multiple scales [3D. 
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